Searching for extremal PPT entangled states 
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We study extremality in various sets of states that have positive partial transposes. One of the 
tools we use for this purpose is the recently formulated criterion allowing to judge if a given state 
is extremal in the set of PPT states. First we investigate qubit-ququart states and show that the 
only candidates for extremal PPT entangled states (PPTES) have ranks of the state and its partial 
transposition (5,5) or (5,6) (equivalently (6,5)). Then, examples of extremal states of (5,5) type 
and the so-called edge states of type (5, 6) are provided. We also make an attempt to explore the 
set of PPT states with ranks (5, 6). Finally, we discuss what are the possible configurations of ranks 
of density matrices and their respective partial transposition in general three-qubit and four-qubit 
symmetric states for which there may exist extremal entangled PPT states. For instance in the first 
case we show that the only possibilities are (4, 4, 4) and (4, 4, 5). 
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I. INTRODUCTION 

Entanglement is the property of states of composite 
quantum systems which is the major resource for quan- 
tum information In the recent years considerable 
interest has been devoted to the entanglement problem, 
i.e., the question of determination whether a given state 
q of a composite systems is entangled or not. This ques- 
tion, although it can be formulated in an elementary 
way accessible to a fresh student of physics, has very 
deep aspects and is related to applications (quantum 
information, quantum communication, quantum metrol- 
ogy) and advanced interdisciplinary fundamental prob- 
lems. The latter concern physics (properties of corre- 
lations in entangled states, preparation, manipulation, 
and detection) , mathematics (theory of positive maps on 
C*-algebras), philosophy (non-locality of quantum me- 
chanics), and computer science (quantum information). 
Operational criteria for entanglement checking exist only 
in very special cases. A famous Peres criterion says that 
every non-entangled (separable) state has a positive par- 
tial transpose (PPT) 0] (see also Ref. [|[). In two qubit 
and qubit-qutrit systems the partial transposition gives 
necessary and sufficient operational criterion of entangle- 
ment jij. In all higher-dimensional systems there exist 
PPT entangled states (PPTES, for the first examples see 
Ref. @). 

Recently, the interest in the entanglement problem has 
evidently become to cease. The reasons for that are 
twofold: on the one hand Gurvits Q has demonstrated 
that the entanglement problem is AP-hard, i.e. finding 
operational criteria in higher dimensions is not likely. On 
the other hand, for low dimensions Doherty et al. [7( 
have formulated a very efficient numerical test employing 
methods of semi-definite programming and optimization. 

Entanglement, quantum correlations and quantum 
measurement theory were in the center of interests of 



the late Krzysztof Wodkiewicz, who made recently re- 
markable contributions to the subject. His works con- 
centrated on the entanglement and correlations in finite 
dimensional and continuous variable systems Q, quan- 
tum non-localit y [9|1 , and implementations in quantum 
optical systems [10]. This paper, in which we reconsider 
the entanglement problem, is dedicated to his memory. 

Our approach is not targeted at operational criteria, 
rather we attempt to characterize and parameterize the 
whole set of extremal PPTES, in particular for 2(g) 4 sys- 
tems. Parameterizing this set allows one to construct 
and describe the set of necessary and sufficient entan- 
glement witnesses 



11 



that fully characterize and de- 
tect all PPTES. For this purpose we need to identify all 
the extremal entangled states in the set of PPT states. 
In earlier works some of the authors have characterized 
other class of PPT entangled states closely related to 
the extremal ones. These are the so-called edge states 
1 3. [la | (for some examples of the edge states see e.g. 
5l Il4l. Il5l fill fl7l|). Quite obviously, edge states are the 
only candidates for the extremal PPTES as every ex- 
tremal PPT entangled state is also of the edge type. In 
the search for these extremal states we arrived at a sim- 
ple criterion which reduces the problem to the existence 
of solutions of a system of linear equations. Later [37| 
we learned that the same result has already been worked 
out in Ref. ^ (see also Ref. 0). Also, in Ref. [H 
the existence of extremal entangled PPTES of bi-rank 
(5, 6) (see below for the explanation of the notation) was 
confirmed numerically. In what follows we provide our 
formulation of the criterion but one of the main pur- 
poses of the paper is to study its applicability in various 
systems as qubit-ququart, three-qubit, and four-qubit 
symmetric density matrices. We also make a step to- 
wards understanding the structure of extremal states in 
2 ® 4 systems. 

It should be noticed that some effort has been recently 
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devoted to a search for extremal PPT entangled states in 
qutrit-qutrit systems. In Ref. [20| 3 ® 3 extremal PPT 
entangled states of type (4_4) were found. Then, the edge 
states provided in Refs. [16|, [1711 o f type (5, 5) and (6, 6) 
were shown to be extremal in [21| and [19| , respectively. 

The paper is organized as follows. In Section |TT] we 
provide our formulation of the criterion. First, we dis- 
cuss a simple observation concerning extremal points of 
an intersection of two convex sets. The observation is 
then applied to convex sets of PPTES and the opera- 
tional extremality criterion is formulated. These ideas 
are then applied in Sec. IIIII to the case of 2 ® 4 sys- 
tems. We discuss general properties of PPTES in such 
systems and then focus our attention on states that have 
ranks of the state and its partial transposition (5, 5) or 
(5, 6) (equivalently (6, 5)), and which are the only candi- 
dates for extremal PPTES. In the case of bi-rank (5, 5) 
we prove that every edge state is extremal in the set of 
PPT states and provide examples of extremal states of 
the type (5,5). Also, examples of (5,6) edge states are 
provided. In Sec. IIVI we make an attempt to explore 
the structure of PPTES of type (5,6). Finally, in Sec. 
fVl we discuss the applicability of the criterion in multi- 
partite systems as three-qubit and four-qubit symmetric 
states. We show also that any three-qubit edge state of 
type (4, 4, 4) is extremal and thus prove extremality of 
the bound entangled state provided in Ref. [l4j. 

Our results are primarily of fundamental interest. 
They shed light on geometry of quantum states in general 
[22| and on the structure and nature of the very complex 
convex set of PPTES, in 2 ® 4 systems in particular. 

II. EXTREMALITY CRITERION 

Here we present our formulation of the criterion (l8j 
for judging if a given element of the set of PPTES is 
extremal. For further benefits it is desirable to set first 
the notation and explain in more details notions which 
have already appeared in the Introduction. 

A. Preliminaries 

Let then g be a bipartite state acting on a product 
Hilbert space TC = C dl ® C d2 . All such states constitute 
a convex set which we denote by T>d 1 ,d 2 - It consists of 
two disjoint subsets of separable and entangled states. 
The distinction between both sets is due to Werner [23| . 
Following [23j, we call separable any state acting on TL 
that can be written as a convex combination of product 
states, that is 

Q = J2Pi$ ) ^e®' p*^ ' &* = 1 > CO 

i i 

with acting on C dj (j = 1,2). If g does not admit 
the above form we say that g is entangled. One of the 



most famous and important tests of separability is based 
on the notion of transposition. Namely, it was noticed 
by Peres that application of the map T®I (hereafter 
called the partial transposition) , where I stands for the 
identity map, to any separable state gives other separa- 
ble state. Simultaneously, when applied to an entangled 
state g the partial transposition can give a matrix (here- 
after denoted by g FA = (T ® I)(g) or g Ts = (I ® T)(g)) 
that is no longer positive. Thus, the transposition map 
may serve as a good entanglement detector. In fact, it 
was shown in Ref. Q that in qubit-qubit and qubit- 
qutrit systems it detects all the entangled states. Inter- 
estingly, there exist states which do not admit the form 
([T]), but nevertheless their partial transposition is posi- 
tive (for the first examples see Ref. H). Entangled states 
that have positive partial transpose cannot be distilled 
with local operations and classical communication to a 
pure maximally entangled states, and therefore are also 
called bound entangled (24[. 

Generally the partial transposition maps the set of 
density matrices T>d 1: d 2 onto other set (containing also 
non-positive matrices) 2? r di> which, as one easily veri- 
fies, is also convex. The intersection of both sets T>^J — 

T^d\,di ^^dj d 2 ' which of course is convex, contains all the 
states with positive partial transpositions. Here we arrive 
at an important notion in the nomenclature connected 
to the set of PPT states. Namely, certain of the PPTES 
that lie on the boundary of 2? PP J are called edge states 
(l3T |. The formal definition says that given PPT entan- 
gled state g € T>d\,d^ 1S an e< ^9 e state if for all e > and 
separable pure states |e, /) = |e)|/) <G C dl ® C d2 it holds 
that g — e|e,/)(e,/| is not positive or its partial trans- 
position is not positive. From this definition it follows 
immediately that any edge state is an entangled state 
that does not contain in its range such a separable vector 
|e, /) that |e*,/) (|e*) denotes the complex conjugation 
of |e)) belongs to the range of its partial transposition. 

We denote further by R(X), K{X), and r{X) the 
range, kernel and the dimension of R(X) (that is the 
rank) of the matrix X, respectively. Moreover, as the 
pair of numbers (r(g), r(g TA )) occurred to be useful in 
a classification of the edge states [l3[ and then com- 
monly utilized in the study of PPT states (see e.g. Refs. 
QJ, EH, [H, QJ, El), hereafter we shall call it bi-rank of 
g. 

B. Extremal points of an intersection 

Let us now discuss what are the extremal points of 
an intersection of two nonempty convex sets <Si and S2 
provided that Si PI S2 is nonempty. Recall that a given 
element x G S is an extremal element of S if it can- 
not be written as a convex combination of elements from 
S which are different from x. There are two classes of 
extremal elements in S\ D S2 - this statement we can for- 
malize as the following (for an illustrative example see 
Fig.©. 
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FIG. 1: The dodecahedron on the left and the known cube 
on the right. A polyhedron consists of polygonal faces, their 
sides are known as edges, and the corners as vertices. Edges 
and vertices are also examples of faces of polyhedron. 



Observation. Let Si and S2 be two convex sets. Then 
the set of extremal points of Si n S2 consists of: (i) the 
extremal points of S± and S2 belonging to Si n S%, (ii) 
possible new extremal points of Si D S2 which belong to 
the intersection of boundaries of Si and S2 and are not 
extremal points of Si and 52. 

To proceed in a more detailed way let us introduce 
some definitions. Let S be some convex set and let x, y 
(x 7^ y) be two elements from S. Then any subset of S of 
the form [x, y] = {Xx + (1 - X)y | x, y £ S, < A < 1} is 
called closed line segment between x and y. Further, let F 
be some subset of the convex set S. Then F is called face 
of S if it is convex and the following implication holds: 
if for any line segment [x, y] C S such that its interior 
point belongs to F then [x, y] C F. Equivalently, F C S 
is a face of S if from the fact x,y £ S and (x + y)/2 £ F 
it follows that x,y £ F. Thus, x £ S is an extremal point 
of S if the one-element set {x} is a face of S. 

This definition formalizes the notion of a face of a con- 
vex polygon or a convex polytope and generalizes it to an 
arbitrary convex set. In this sense, vertices of a polytope 
are zero-dimensional faces (extremal points), edges are 
one-dimensional faces and the traditionally understood 
faces are two-dimensional. The interior points form a 
two-dimensional face (see Fig. [I}. On the other hand, 
any point on the boundary of a closed unit ball in R™, 
i.e., any point of the unit (n — 1) -dimensional sphere, is 
a face (and an extreme point) of the ball. 

Let us now discuss what are the main conclusions fol- 
lowing the above observation. Namely, if some convex 
sets Si and S2 intersect in a point x along an affine man- 
ifold of non-zero dimension then x cannot be extremal. 
Therefore x belongs to the second class of extremal points 
in Si H S2 if and only if the intersection of the relevant 
faces is exactly x, i.e., the faces are transversal at x (see 
Fig. [2]). 

Let us also point out how the above observation can be 
utilized to prove that a given x is not an extremal element 
of Si n 52- Namely, if x belongs to the intersection of 
some faces of Si and S2, we can consider other element 
x±ey with e denoting some arbitrarily small positive real 
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FIG. 2: A simple example visualizing the observation formu- 
lated in Sec. Ill Bl The se of extremal points of the rectangle 
A'ECF resulting from an intersection of rectangles ABCD 
and A'B'C'D' consist of some of the extremal points of both 
rectangles or the ones that appear as intersections of their 
respective edges. More precisely, the points A' and C are 
extremal points of the rectangles A'B'C'D and ABCD, re- 
spectively. The Points E and F result from intersection of 
the edges of these rectangles, namely, A'B' and BC, and CD 
and D'A', respectively. Both the intersections give manifolds 
consisting of a single point (the respective edges are transver- 
sal). 



number and y is an element of the affine space in which Si 
and S2 are convex subsets. If such element belongs to the 
intersection of the mentioned faces the intersection has a 
non-zero dimension and x cannot be an extremal element 
of Si<lS2- The key point of this remark is that to generate 
this manifold we do not have to consider elements of the 
convex subsets Si and S^- This, as we will see in the 
next subsection, can be applied to convex sets of density 
matrices with positive partial transposition. 



C. Criterion for extremality in convex sets of 
quantum states with positive partial transposition 

In the context of T>^J the observation of Sec. Ill Bl 
says that the set of its extremal points is the sum of two 
disjoint subsets. The first one consists of pure extremal 
states of £>di,d 2 which are simultaneously extremal points 
of T> v di d . The second subset consists of the extremal 
points that could appear from intersection of boundaries 
of ,d 2 an d d 2 anc ^ are n0 ^ extremal points of both 
these sets. If a given g is neither an extremal point of 
^di,da nor °f f di d 2 then we have to check (assuming 
obviously that it is an edge state) if the faces to which 
g and g TA belong are transversal. It can be seen that 
this reduces to the problem of solving a system of lin- 
ear equations. For this purpose we utilize the technique 
mentioned in the preceding section of adding or subtract- 
ing small elements not necessarily belonging to 2?d 1 ,d 2 
or Z>£ d2 . As we will see below for this aim we can 
use Hermitian matrices H such that R(H) C R(g) and 
R(H rA ) C R(g rA ). To proceed more formally let us no- 
tice that, as we know from Ref. (25|, the face of 2?^ <j a at 
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a given g (which is not of full rank as only such density 
matrices constitute faces of T> dltd2 ) consists of matrices 
of dimensions d\d 2 x d\d 2 with kernels coinciding with 
the kernel of g (this holds also for partial transpositions 
of density matrices provided they are positive). On the 
other hand, all Hcrmitian matrices of rank k constitute k 2 
dimensional subspace of the (did2) 2_ dimensional Hilbert 
space V (equipped with e.g. the Hilbert-Schmidt scalar 
product) of all d\d% x d\d% Hermitian matrices. To check 
if the intersection V\ fl V 2 of two subspaces V\ and V 2 of 
dimension k 2 and I is nonempty, one has to solve the 
system of 2(did2) 2 — k 2 — I 2 linear equations for (ditfe) 2 
variables. A nontrivial, nonzero-dimensional, solutions 
exist always if (d\d2) 2 < k 2 + I 2 (in the case of density 
matrices, due to the normalization, we add one to the 
left-hand side of this inequality) . If our density matrix 
g has a bi-rank (k, I) admitting the existence of such so- 
lutions, it cannot be extremal in T> d J 2 - one can always 
find a Hermitian matrix H £ Vi fl V 2 with R{H) C R(g) 
and R(H Fa ) C R(g TA ) such that g(e) = g±eH is positive 
and PPT. What is important, for small enough e > the 
matrix g(e) has the same bi-rank as the initial state g. 
From the fact that the face of T>d 1: d 2 at a given g (which 
is not of full rank as only such density matrices consti- 
tute faces of 2?di,d 2 ) consists of matrices of dimensions 
d\d 2 x d\d 2 with kernels coinciding with the kernel of g 
it follows that g(e) has to belong to the same intersection 
as g and in such case g cannot be extremal in T> FF J 2 . In 
this way we finally arrive at the criterion. 

Criterion. Let g £ T>d ljd2 be PPTES with the bi-rank 
(k,l). If there exists a non-trivial, i.e. non-proportional 
to g, solution of the system of linear equations described 
above, i.e. equations determining Hermitian matrices H 
with R(H) C R(g) and R(H Fa ) C R(g TA )), then g is not 
extremal in T> FF J . 

It should be emphasized that the above technique of 
adding and subtracting Hermitian matrices and the cri- 
terion following it were already provided in Ref. :18|. 

In conclusion, to prove that some g is not extremal in 
the set of PPT states one can always try to "generate" the 
nonzero-dimensional manifold by adding or subtracting 
eH with H being a general Hermitian matrix H and e 
some sufficiently small positive number. This, however, 
has to be done in such a way that the resulting state 
f?( e ) = g±eH 7 as well as its partial transposition have the 
same kernels as, respectively, g and g Fjl . More precisely 
we need to have R(H) C R(g) and R(H Fa ) C R(g FA )). 
In this way we remain on the same faces of T>d t ,d 2 and 
T>j d2 (as the partial transposition of PPT state is still 
a legitimate state then the observation above works also 
for partial transpositions of density matrices) as g itself. 



III. THE SET OF 2 g> 4 STATES WITH THE 
POSITIVE PARTIAL TRANSPOSITION 

Here, we sketch shortly the overall picture of qubit- 
ququart states with respect to the extremality in the set 



of PPT states. First we recall literature results from 
which it follows that the only possible cases of PPT en- 
tangled states that could be extremal are the ones with 
bi-rank (5,5) and (5,6) (and equivalently (6,5)), and 
(6,6). Next, basing on the observation from Sec. HH we 
will exclude the case of bi-rank (6,6). Then, we show 
that any edge state of bi-rank (5, 5) is extremal in T> FPT 
which implies that the famous examples of 2 <g> 4 PPT en- 
tangled states provided in Ref. [f| are extremal. Finally, 
we present a two-parameter class of edge states that de- 
pending on the parameters have bi-rank (5,6) or (5,5), 
giving at the same time next example of extremal (5, 5) 
states. 



A. The general situation 

Let us here discuss shortly why the only possible cases 
in which there could be PPT entangled extremal states 
are these of the bi-rank (5, 5), (5, 6) (equivalently (6, 5)), 
and (6,6). At the very beginning we can rule out the 
cases in which either g or g FA are of full rank since then 
the density matrix lie in the interior of the set of PPT 
states (see e.g. Ref. [25j). 

To deal with the remaining cases we can assume that 
the given g is supported on full C 2 (8) C 4 . Otherwise we 
would deal with a state that is effectively qubit-qutrit 
state and one knows that in this case there are no PPT 
entangled states Q ■ On the same footing we can assume 
that there is no separable vector |e, /) in the kernel of g. 
Otherwise, it would mean, according to Ref. .12], that 
g could be written as a mixture of a projector onto the 
separable vector \e ,h) (with (e |e) = 0, and a certain 
\h)) and some PPT state supported on C 2 <£> C 3 , and 
thus separable. In conclusion, in such a case g could be 
written as a mixture of separable states and thus would 
be separable itself. 

We know from Ref. [12] that if g is supported on C 2 ® 
C N and is entangled then r(g) > N. This means that 
in what follows we can assume that r (g), r(g FA ) > 4 as 
otherwise we would deal with a separable state. Then, 
following Ref. [II] , we divide our considerations into two 
cases, 

(i) r{g)+r{g FA ) < 12, 

(ii) r(g)+r(g FA ) > 12. 



In the second case it was shown in Ref. 12] that there 
always exists a separable vector |e, /) in the range of g 
for which |e*,/) € R(g FA ). Consequently, we may sub- 
tract this vector from g obtaining another (unnormal- 
ized) state 

g = g-v\ej)(ej\, (2) 

where r) = min{r/o,?7o} with Vo an d Vo gi ven by t]q — 
{{e,f\g-_ X \e,f))- X and % = {{e% f\{g FA )^\e* , f))~\ 
Here g 1 denotes the so-called pseudoinverse of g, i.e., 
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inverse of g only on the projectors corresponding to its 
nonzero eigenvalues. 

Normalizing g we infer from Eq. ([2]) that g can be 
written as 



tremal in T>\ 4 T . 



This means that it admits the form 

PPT 



Q= (1 - V)Q + V\e,f)(e,f\, 



(3) 



with < r\ < 1 (this is because both parameters rjo and 
rj satisfy < r] ,rj < 1). The important fact here is 
that the above procedure of subtraction product vectors 
preserves positivity of partial transposition of g. More 
precisely, if g r - A > then also g TA > 0. All these mean 
that if r(g) + r(g rA ) > 12 we are always able to write 
g as a convex combination of two states with positive 
partial transposition. Thus in this case there are no PPT 
entangled extremal states in the convex set T> f 4 T (we see 
that even edge states do not exist in these cases). 

So far we have ruled out most of the cases with re- 
spect to the bi-rank. What remains are the states satis- 
fying r(g) > 4, r(g FA ) > 4, and r(g) + r(g FA ) < 12. In 
the cases of (5,7) and (7,5) it was shown in Ref. [26[ 
that there exists a product vector |e, /) E R{g) such 
that \e*,f) € R(g rA ). This, according to what was said 
above, means that also among states with bi-rank (5, 7) 
or (7, 5) one cannot find PPT entangled extremal states. 
The case of bi-rank (6, 6) can be ruled out by the cri- 
terion. Here k = I = 6 and therefore one has a system 
of 56 linear homogenous equations for 63 variables (one 
is subtracted due to normalization). According to what 
was said in Sec. Ill Cl this means that all (6, 6) PPTES are 
not extremal as for any such state there exist a Hermi- 
tian H satisfying R(H) C R(g) and R(H Fa ) C R(g FA ). 
Thus an intersection of the corresponding faces of 2?2,4 
and T) T 4 is thus a nonzero-dimensional manifold. 



B. The case of bi— rank (5, 5) 

The remaining cases in which we can look for the ex- 
tremal PPT states are (5,5) and (5,6). In the first case 
of the bi-rank (5, 5) examples of extremal PPT entangled 
states can be provided. For instance, we can show that 
the states found in Ref. [5j, 
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are extremal for 6 € (0, 1) (for 6 = 0, 1 the state is sepa- 
rable). For this aim let us prove the following theorem. 

Theorem. Any 2 4 edge state with bi-rank (5, 5) is 
extremal in the set of PPT states T>\ 4 T . 

Proof, (a. a.) Let us consider an 2 ® 4 edge state a 
with the bi-rank (5, 5) and assume that it is not ex- 



(7 = \o\ + (1 — A)(72 with some A G (0, 1) and at € 2}ij 4 
(i = 1,2). Since a is edge state it is clear that both 
the density matrices Oi have to be entangled. Also, it 
is easy to see that the conditions R(oi) C R{a) and 
R{o-f A ) C R(a TA ) hold for i = 1,2. As a conclusion 
one of these density matrices, say ci, can be subtracted 
from a with some small portion e > in such a way 
that the resulting matrix, as well as its partial transpo- 
sition, remain positive. More precisely, we can consider 
the state a = (a — e<Ji)/{\ — e) which obviously has the 
rank r(a) < 5. The parameter e can be set so that either 
r(a) = 4 or r(g r ) = 4. Consequently, by virtue of re- 
sults of Ref. [12| where it was shown that any PPT state 
supported on 2 (g) 4 with rank 4 is separable, the latter 
means that a is separable. As a result the initial state a 
can be written in the form 



(1 



ejcr + eoi, 



(5) 



which means that a is a convex combination of a sepa- 
rable and an entangled state. This, however, is in con- 
tradiction with the assumption that a in an edge state. 
Thus a has to be extremal what finishes the proof. | 

Notice that the above statement remains true if one 
relaxes the condition of being an edge state to being a 
PPT entangled state. That is, all PPT entangled states 
of bi-rank (5, 5) are extremal. 

It was shown in Ref. [j| that <r(6) for 6 € (0, 1) are edge 
states and therefore, due to the above theorem, they are 
extremal in the set of 2 ® 4 PPT states. 

Let us remark that all the cases satisfying r(g) + 
r(g FA ) > 12, as well as the case of the bi-rank (5,7) 
(and equivalently (7,5)) can be also ruled out using the 
criterion. 



C. The case of bi— rank (5,6) 

The last case we need to deal with is the one with 
the bi-rank (5,6). It was numerically confirmed in Ref. 
[l8l ] that there exist (5, 6) extremal PPT states in T> PPT , 
however, explicit examples are still missing. 

In what follows let us consider the family of 2(g) 4 PPT 
entangled states with r(g) = 5 and r(g TA ) = 6 which are 
edge but not extremal in the set 2?2 4- The family is of 
the form 
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The vector of eigenvalues of g(a, t) reads 
1 



X(g(a,t)) = 



2(1 + a) 2 

x (2a, 1 + a 2 , 1 + a 2 , a(l - t), a(l + t), 0, 0, 0) . (7) 

Thus, for a > and — 1 < t < 1 the rank of g(a,t) is 
exactly five. Let us now look at the partial transposition 
of g, which is given by 



[g(a,t)} 1 
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2(2 + a + a- 1 ) 
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(8) 

The vector of eigenvalues of g TA is of the form 

A ^ a '^^2(iT^ 

x (1 - a, 2a, 2a, 1, a(a - t), a(a + 1), 0, 0) . (9) 

For a > t, a + t > 0, and a < 1 one has [g(a,t)) rA > 0, 
hnece g(a, t) is a PPT state of the bi-rank (5, 6) whenever 
< a < 1 and \t\ < a. For t = a it has the bi-rank (5, 5) 
and is an edge state (one can check it using the method 
presented in Ref. Q). Due to the theorem from Sec. lIII Al 
it is also another example of an extremal state in T>\ 4 T . 

To prove that the g(a, t) is entangled and thus bound 
entangled we may use the already mentioned range cri- 
terion of Ref. After some algebra one finds that all 
the product vectors in the range of g(a, t) can be written 
in the form 

A(a,t)(l; a) ® (a 3 ; -a 2 /a; a/ a; 1) (aeCU {oo}), 

(10) 

with A being some function of a and t. We used here the 
fact that any vector from C 2 can be written as |0) + a|l) 
with a G C U {oo} end employed the convention a|0) + 
6|1) = (a; b) together with a similar one for vectors in C 4 . 

A similar analysis shows that the product vectors in 
the range of g TA are of the form 

(1; a) (g> (B(a, a, t), —aB(a, a, £),C(a, a, t), —aC(a, a, t)) 

(aeCU{oo}), (11) 

with B and C denoting some functions. Simple calcula- 
tions show that one cannot find such a product vector 
|e,/) from R{g(a,t)) that |e*,/) g R(g(a,t) FA ). More 
precisely, one cannot find such functions B and C in Eq. 
(fTTj) that the vector A(a, t)(l; a*) ® (a 3 ; -a 2 /a; a/a; 1) 
can be written in the above form. Thus, for any a and t 
satisfying < a < 1 and \t\ < a, the matrix g(a,t) is a 



PPT entangled state. Moreover, the same argument tell 
us that it is also an edge state. 

Now, one immediately finds that g(a, t) cannot be ex- 
tremal in the set of 2 ® 4 PPT states. This is because 
for a fixed a £ (0, 1), t E (—a, a) can always be expressed 
as t = (1 — A)ii + Xt2 with A e (0, 1) and U € (—a, a) 
(i = 1, 2), which means that for any a £ (0, 1) the state 
g(a, t) can be written as the following convex combina- 
tion 



p(a, t) = (l- X)g(a, ti) + \g(a, i 2 ). 



(12) 



IV. 



EXPLORING PPT STATES OF BI-RANK 

(5,6) 



We continue our considerations in the case of bi-rank 
(5, 6), however, in the most general situation. As the dis- 
cussion presented in this section is a bit more mathemat- 
ically demanding it is helpful to outline first its contents 
before going into details. The general aim of this section 
is to ex plor e PPT states with bi-rank (5,6). It is stated 
in Ref. [lj| that there exist extremal PPTES of this type, 
however, the general structure of this class of states is not 
known and the main purpose is to make a first attempt 
to understand it. For this purpose, using the so-called 
canonical form of 2(g) TV states, we determine all the prod- 
uct vectors in ranges of g and g TA . Then, using the lat- 
ter and following the criterion, we propose for any (5, 6) 
PPT state such a Hermitian matrix H that R(H) C R(g) 
and look if R(H Fa ) C R(g TA ). This brings the problem 
of extremality to the problem of solving a set of linear 
equations. Although we cannot provide a solution to this 
system in the general case, we can provide various sets 
of conditions for B (see Eq. (fTS")) ) provided that the so- 
lution has some particular form (and the corresponding 
state is not extremal). 



A. Separable vectors in ranges of g and q Ta in the 
case of (5, 6) 

Let us consider an arbitrary 2 ® 4 state. It can be 
brought to the so-called canonical form (see e.g. Ref. [26|, 
H3 and Ref. [Hi in the context of 2 ® 2 ® N states). To 
see it explicitly let us notice that any such state has the 
form 



A B 
fit c 



(13) 



where A, B, and C are some 4x4 matrices and due to 
the positivity of g one knows that A and C have to be 
positive as well. Since we assume that r(g) — 5, we can 
also assume that the matrix C is of full rank. Otherwise, 
one sees to see that the vector |1, /), where |/) £ K(C), 
belongs to the kernel of g. In view of what was said before 
this means of course that g would be separable. 
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Thus, assuming that r(C) — 4 we can bring g to its 
canonical form by applying the transformation g i— * (l2<8> 
C~ 1/2 )g(l 2 ® C~ 1/2 ), which finally gives 



= 



A 1 B' 
B't 1 4 



(14) 



where by Id we denoted d x d identity matrix. It should 
be emphasized that the above transformation does not 
change the rank, extremality as well as separability prop- 
erties of a given g. Thus from our point of view both 
states g and (I2 &> C~ 1 ^ 2 )g(t 2 <B> C -1 / 2 ) are completely 
equivalent. 

We can simplify the form (fT4|) even further. The pos- 
itivity of p' implies A 1 = B'B'* + A, where A is some 
positive matrix. Consequently, omitting all the primes 
we can write 



BB^+k B 
St 1 4 



(15) 



where B is some (in general non-Hermitian) matrix. Fur- 
thermore, we can assume that det(i?) = 0, which means 
that r(B) = 3 and det(B T ) = 3 and therefore there exist 
such \4>) and \4>) from C 4 that B\<f>) = and B^\4>) = 0, 
respectively. This can be done by changing the basis in 
the first subsystem. Precisely, we can apply the transfor- 
mation C/(7, 8) (8 I4 with [7(7,(5) given by 



u(rr,8) 



7 8 
-8* 7* 



( 7 ,<SeC\{0}), (16) 



to g. Then it suffices to solve the equation det(B) = 
which always has solutions. Finally, A is some pos- 
itive matrix of rank one and thus can be written as 
A = A 1 |A 1 )(A 1 |. 

Performing the partial transposition with respect to 
the first subsystem we get 



9 TA = 



BBi - 
B 



A St 

1 4 



Et£-t 
B 



A St 
1 4 



(17) 



with A = [B, fit] + A = A 2 |A 2 )(A 2 | + A 3 |A 3 )(A 3 |. This 
form is similar to the one of g with the difference that 
A is a combination of two projectors as we want to have 
r(g FA ) = 6. 

Let us now find the general form of separable vectors in 
the ranges g and g TA and discuss their properties. Due to 
the similarity of the forms of g and g TA (see Eqs. (|T5|) and 
(|17p). it suffices to make the calculations for the range 
of g. As any (unnormalized) vector from C 2 ® C 4 can 
be always written as \0,g) + \l,h), this can be done by 
solving the system of equations g(\0,g) + \l,h)) = \e)\f) 
with |e) £ C 2 and \g), \h), and |/) denoting some vectors 
from C 4 . As previously we can utilize the fact that any 
vector |e) from C 2 can always be written as |e) = |e(a)) = 
|0) + a|l) with a £ C U {00} (for the sake of simplicity 
we forget about normalization, however we keep in mind 



that |e(0)) = |0) and |e(oo)) = |1)). The resulting set of 
equations reads 



(BBt+A)\g)+B\h) = \f), 
Bi\g) + \h) = a\f). 



(18) 



After some algebra one finds that the solution is of the 
form 



|/) = a|/(a))=a(l 4 -aS)- 1 |A 1 ) ) 



(19) 



with a = Ai (Ails')). Except for the solutions of det(l4 — 
aB) — 0, the inverse of I4 — aB exists for all a £ C. 
The same reasoning in the case of g TA gives 



\f(a,b,c)) = (l-aB^- 1 (b\X 2 )+c\X 3 )). 



(20) 



with appropriate constants b and c. Therefore the sep- 
arable vectors from R(g) and R(g TA ) are of the form 
a\e(a))\f(a)) and \e(a))\f(a,b,c)), respectively. With 
some additional effort we can determine a more detailed 
form of the separable vectors in the range of g proving 
that the general form of |/(a)) is 



!/(")> 



1 



W{a) 



(W 1 (a);W 2 {a);W 3 (a);W i {a)), (21) 



with W(a) and Wi(a) (i — 1,...,4) being some poly- 
nomials in a of degree at most three. To see this ex- 
plicitly we notice that since r(g) = 5 there exist three 
linearly independent vectors \4>i) in K(g) orthogonal to 
|e(a))|/(a)). Writing them in the form 



|<MH0)|$ 0) ) + |i)lO, (22) 

we get from the orthogonality conditions 

(^ 0) |/(«)> + a<^ 1) l/(«)> = 0, (i = l,2,3). (23) 

By decomposing \f(a)) — Y^k=i /fe( a )l a fc) m some basis 
{|afc)} we get the following system of three homogenous 
equations 



$>(«)(<c>*>+« 



r(i) 1 



o fc ) =0. 



(24) 



We can always fix the value of one of e.g. by putting 
/4(a) = 1, and obtaining 



fc=i 



m 1 



^/ fc (a)((^ 0) |a0+«(O0j = -WW 

-all 11 !"!), 



(25) 

for i = 1,2,3. Now we deal with a system of three 
inhomogenous equations (the inhomogeneity is nonzero 
as one can always find such basis in C 4 that at least 
one of the scalar products on the right-hand side of 
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(|25| is nonzero). Solving the system we get the pos- 
tulated form ([2"T]) . Using a little bit more sophisticated 
reasoning one may also prove that the polynomials Wi 
(i = 1,2,3,4) are linearly independent and therefore at 
least one of them must be of the third degree. This 
also means that all the separable vectors in R(g) can be 
brought to the form (l;a) ® (1; a, a 2 ; a 3 ). Indeed, since 
Wi (i = 1,2,3,4) are linearly independent there exist 
matrix V with det(U) ^ such that 1 2 ® Vgt 2 ® has 
all the separable vectors in its range of the above form. 



B. Density matrices of bi— rank (5,6) that are not 
extremal 



Knowing the general form of product vectors from 
R(g) and R(g TA ) of a general qubit-ququart state we 
can now provide some conditions under which the (5, 6) 
PPT states are not extremal. For this purpose let us 
consider the following Hcrmitian operator 

H = oia 2 |e(ai)Xe(a 2 )| <8 |/(ai))</(a 2 )| + h.c. (26) 

with a\ ^ a 2 . Its partial transposition with respect to 
the first subsystem reads 

H Ta = oiOa|e(a;)Xe(a£)| (8 |/(a0)</(a 2 )| + h.c. (27) 

One sees that R(H) C R(g), still, however, we need to 
check if R(H Ta ) C R(g rA ). The assumption ot\ ^ a 2 is 
important since for ct\ — a 2 one could subtract a product 
vector from g without spoiling its relevant properties, 
hence g would not be an edge state. 

We need to find such a±, a%, a 2 , and a 2 
that ai|e(a 2 ))|/(o!i)) = Je(a 2 ))|/(a|, ai, &i)) and 
a 2 \e(al))\f(a 2 )) = \e(al))\f(al, b 2 , c 2 )) as both these 
vectors belong to the range of g TA . These conditions 
simplify to ai|/(ai)) = \f(a 2 , h, Ci)) and a 2 |/(a 2 )) = 
\f{a\, b 2 , c 2 )), respectively which, in turn, can be brought 
to a system of linear equations 

(J-a x B)|a;)=oi|Ai), (28) 

(I-a* 2 B)\x)=b 1 \X 2 )+c 1 \X 3 ), (29) 

where \x) is a vector from C that has to be determined. 
We can introduce three vectors \fi) (i = 1, 2, 3) orthogo- 
nal to | Ai) and two vectors \vi) (i = 1,2) orthogonal to 
|A 2 ) and |A 3 ). Projection of Eqs. (J2HJ) and onto the 
vectors |/j) (i = 1,2,3) and |A X ) leads to the following 
systems of equations 



(30) 



f (fi\I-aiB\x) = 0, i = 1,2,3, 
| (Ai| J - aiB\x) = ai, 

and onto the vectors \vi) (i — 1,2) and |A 2 ), and | A3) to 
(vi\I-a* 2 Bi\x) =0, i = 1,2, 
(\ 2 \I-a* 2 Bt\x)=b 1 + ci(\ 2 \\ 3 ), (31) 
I (A 3 |/-a5Bt| a; ) = 6 1 (A3|A 2 )+c 1 . 



The same procedure applied to the equation a 2 |/(a 2 )) = 
\f{al, b 2 , c 2 )) leads to an analogous system of equations 

f (fi\I-a 2 B\y)=0, i = 1,2,3, 
\ (Ai|7 - a 2 B\y) = a 2 , 



(32) 



{Vi\l-a^\y)=0, i = 1,2, 

(A 2 |/-aJSt|y)=& 2 + C2 (A 2 |A 3 ), (33) 

(A 3 |J-aIBt| y ) =6 2 (A 3 |A 2 ) + c 2 . 

Together, the homogenous equations in (|3H)) . PT|) . (|32p . 
and (|33p give ten equations for eight unknowns variables 
Xi and j/i (i = 1 . . . , 4) being coordinates of \x) and \y) 
in some basis in C 4 . From the remaining inhomogenous 
equations one could determine a,, bi, and Cj (i = 1,2). 
The most natural and, in our opinion, easiest way to 
solve all the homogenous equations is to consider them 
separately for |x) and \y). For this aim let us rewrite the 
homogenous equations of and (f3"Tj) as 



{fi\I-a x B\x)=Q, i = 1,2,3, 
(vi\I - a%Bi\x) =0, i = 1,2. 



(34) 



They constitute five linear equations for |a;), which im- 
plies that two determinants have to vanish; we take them 
to involve the first three equations, together with one 
of the remaining two. In effect each determinant has a 
form of a third order polynomial in ct\ and first order 
in a 2 . Such a system of two polynomials has several 
solutions for the pairs (a\,a 2 ) with a% 7^ a 2 (if there 
were a solution a% = a 2 , we could subtract the projec- 
tor on the corresponding product vector from the state 
g and the partially transposed one from the partially 
transposed state keeping the positivity of both, which 
is incompatible with g being an edge state) . The remain- 
ing equations, (fi\I — a 2 B\y) = 0, with i = 1,2,3 and 
(vi\I — a\B^\y) = 0, with i = 1,2 we obtain by replacing 
in d3H) \x) — > \y), and ct\ — > a 2 . Concluding, there al- 
ways exist such pairs of ct\ and a 2 that the set (|34l) has 
a nontrivial solution. Now, if the analogous set of equa- 
tions for the vector \y) has a solution for the particular 
ct\ and a 2 solving (|3~4"|) , but interchanged a\ <-* a 2 , then 
we have a solution to the initial problem. We cannot find 
the general solution to this problem and due to results of 
Ref. [l8| it is even impossible. Thus, we consider some 
particular cases with respect to a% and a 2 , and assum- 
ing that the solution is of some particular form we derive 
certain conditions that have to be satisfied by the matrix 
B (see Eq. ([IS])). 

For ai = 00 and a 2 = the homogenous equations 
from PCJ ]) . (|5T] ) . ([52 ]) , and ([55 ]) reduce to 



and 



(/i|B|a;)=0, i = 1,2,3, 

(Vi\x)=0, i = l,2, 

{fi\y) = 0, i = l,2,'S 

(» i |flt| y )=o, i = l,2. 



(35) 



(36) 
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From the first set we see that \x) = ^|A2) + 77IA3), 
while from the second one that \y) = /J-\Xi). So in the 
case of ol\ = 00 and a 2 = we are able to find |a;) 
and \y), however the matrix B has to satisfy the con- 
ditions (/i|B(£|A 2 ) +J?|A 3 )) = with i = 1,2,3 and 
(wi|St|Ai) = for i = 1,2. 
In the case a\ = and ^ a 2 ^ 00 we have 



(fi\x) = 0, 



1,2,3, 



and 



{v t \I - a* 2 B^\x) = 0, i = 1,2, 



(fi\I-a 2 B\y)=0, i = l,2,3, 



i\v) = o, 



i = 1,2. 



(37) 



(38) 



Here, analogously to the previous case, one has \x) = 
A*|Ai) and \y) = £|A 2 ) —I— 77 1 A 3 ) . Inserting the above forms 
of \x) and \y) to Eqs. (|3T)l and ([55)) we obtain five 
conditions for the matrix B. More precisely, we have 
(wi|J-a|St|Ai) = withz = 1,2 and (/i|J-a 2 S(^A 2 ) + 
^A 3 ))=0 fori = 1,2,3. 

Finally in the case of a,\ = 00 and ^ a 2 7^ 00 the 
systems of equations reads 



(/ i |B|aj)=0, 



1,2,3, 



and 



(Vi\I - a* 2 B^\x) =0, 



(/ i |/-a a B|tf)=0, 
(WilBtly) =0, 



* = 1,2, 



1,2,3, 
1,2. 



(39) 



(40) 



All the provided equations in each of the above cases are 
obviously equipped with conditions of vanishing determi- 
nants. 



V. MULTI-QUBIT PPT EXTREMAL STATES 

It seems interesting to extend the criterion to the case 
of many-qubits states. In what follows we discuss how it 
can be applied to the general three and symmetric four 
qubit states with all partial transposes positive (hereafter 
PPT states). As we will see in both cases we deal with 
an intersection of three convex sets. Let us then discuss 
how extremal points appear when three convex sets Si, 
S2, and S3 intersect provided that Si n 5 2 H S3 is not 
an empty set. Following the observation formulated in 
the previous sections (we can apply it recursively to the 
pair of sets S% and S2 , and then to the pair Si f~l S2, and 
S3) one sees that in general the set of extremal points of 
Si n S2 n S3 consists of (i) extremal points of one of these 
three sets, (ii) points that appear as an intersection of 
faces of two of these three sets (but are elements of the 
interior of the remaining set) and are not extremal in 
any of them, (iii) points that appear as an intersection of 
faces of all the sets. 



A. Three qubits 

At the beginning let us notice that separability prop- 
erties in three-qubit systems were thoroughly studied in 
Ref. (28|. Then all these states were classified in Ref. 

Let now T>2 2 2 be the set (obviously convex) of all 
three-qubit states having all partial transposes positive. 
It is clear that positivity of all partial transposes is equiv- 
alent to positivity of arbitrary two elementary subsys- 
tems, say A and B (positivity of the remaining partial 
transposes g r ° , q Tac , and q Tbc is then guaranteed via 
the positivity of the global transpose). Consequently, 
O222 is intersection of three convex sets, namely, the set 
of all three-qubit density matrices denoted by 2? 2 ,2,2 and 
two sets T>2 2 2 an d f 2 2 2 containing partially transposed 
three-qubit density matrices with respect to the A and 
B subsystem, respectively. 

By k, I, and m we denote ranks of g, q Ta , and q Tb , 
respectively (generically 1 < k, I, m < 8, however, we 
can exclude the rank one cases as they are already ex- 
tremal). First, let us recall that it was shown in Ref. 
[28[ that every three-qubit state g with r(g) < 3 sup- 
ported on C 2 <g> C 2 <g> C 2 is fully separable. Thus we can 
assume that r(g), t{q Ta ), r(g Tn ) > 4 as the same argu- 
ments can be applied to q Ta and q Tb . Then, it should 
be noticed that any three-qubit state can be treated as 
2 eg) 4 state with respect to the partitions A\BC, B\AC, 
and C\AB. This, due to results of Ref. [l2j means that 
every three-qubit PPT state with rank r(g) = 4 is sepa- 
rable across all the above cuts. However, an example of 
a state separable across A\BC, B\AC, and C\AB cuts 
but nevertheless entangled was provided in Ref. [3] (the 
bound entanglement constructed from the so-called un- 
extendible product basis), which as we will see later is 
extremal in T>\ 2 ^- Thus in what follows we can put 
k,l,m> 4. 

Now, we can follow analogous reasoning as in the bi- 
partite case and obtain the inequality fc 2 + 1 2 + m 2 > 65, 
which if satisfied implies existence of such a Hermitian 
matrix H that g(e) = q ± eH is a density matrix for 
some nonzero e with all partial transposes positive and 
the same ranks k,l,m (and thus g of such ranks cannot be 
extremal) . Below we consider only such triples that sat- 
isfy k < I < m as the remaining ones one obtains by per- 
muting the ranks fc, I, m. It follows immediately from this 
inequality that there do not exist PPT extremal states 
of the second kind (see the discussion at the beginning of 
Sec. rV)) . This is because if one of the ranks is maximal 
the inequality becomes i 2 + j 2 > 1 = fc, /, m) and is 
always satisfied. On the other hand, application of this 
inequality to the cases when all ranks are not maximal 
allows to exclude the following cases (4,4,6), (4,4,7), 
(4,5,5), (4,5,6), (4,5,7), (4,6,6), (4,6,7), and (4,7,7). 
As a result the only cases that remain are (4, 4, 4) and 
(4, 4, 5). In the first case an example of extremal PPTES 
was given in [LH . To see it explicitly one may prove a sim- 
ilar statement to the theorem from Sec. IIV Al Namely, 
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adapting its proof and utilizing the results of Ref. [28] 
we get the following theorem. 

Theorem 2. Any three-qubit PPT entangled state sup- 
ported on C 2 ® C 2 ® C 2 with ranks r(g) = r(g rA ) = 
r(g TB ) = 4 is extremal in X>2 2 2- 

Since, by construction the three-qubit bound entangled 
state from Ref. [l4[ is an edge state it has to be extremal 
in V\H of type (4,4,4). 

B. Four qubits in a symmetric state 

In the case of four qubits we restrict our discussion only 
to the states acting on 6>((C 2 )® 4 ), which is a subspace 
of (C 2 )® 4 consisting of all pure states symmetric under 
permutation of any subset of parties (let us call such 
states symmetric). General ./V-qubit symmetric states 
were recently investigated in numerous papers (see e.g. 
Refs. HSIMIIIIII). 

One of the motivations to study four-qubit symmetric 
states comes from a long-standing open question whether 
there exist such states with positive partial transpose 
with respect to all subsystems. If it were so there would 
also exist extremal PPT symmetric states. Notice that 
such states, however, with higher number of qubits were 
found very recently in Ref. [331 ] . In what follows we would 
like to discuss what are the possible cases with respect 
to ranks of respective partial transposes for which there 
could exist entangled PPT extremal symmetric states. 

Notice relaxing the PPT condition to two-party trans- 
poses only as for instance g TAB > and leaving the possi- 
bility that single-particle partial transposes are nonpos- 
itive, it is possible to find four-partite bound entangled 
symmetric state [34| which even allows for maximal vio- 
lation of Bell inequality [35] . Interestingly, this state was 
very recently realized experimentally [361 ]. 

Let us discuss firstly some of the properties of mixed 
states acting on iS((C 2 )® 4 ). Clearly, it is spanned by 
the vectors |0)® 4 , |1)® 4 , \W), \W), and (l/\/6)(|0011) + 
|0101) + 1 1001) + |1010) + |1100) + |0110)), where \W) 
is the so-called W statejW) = (1/2)(|0001) + |0010) + 
1 0100) + 1 1000)) and \W) is obtained from \W) by re- 
placing zeros with ones and vice versa. Consequently, 
all symmetric four-qubit density matrices have ranks at 
most five. It also follows from the above that there are 
only two relevant partial transposes defining set of PPT 
symmetric states. Namely, positivity of partial trans- 
pose with respect to a single subsystem (two subsystems) 
means positivity of partial transposes with respect to all 
single subsystems A, B, C, D (arbitrary pairs of subsys- 
tems AB, AC, etc.). The three-party partial transposes 
are equivalent to the single-party partial transposes via 
the global transposition. Hence, as in the case of thrcc- 
qubits, the set of PPT symmetric states is an intersection 
of three sets, namely, the set of symmetric density ma- 
trices and the sets of partial transpositions of the latter 
with respect to the A and AB subsystem, respectively. 
So, there are three relevant parameters, r(g), r(g rA ), and 



r(g TAB ). Using the fact that we deal with the symmetric 
states we can easily impose bounds on these three ranks. 
We already know that r{g) < 5. Then, since S((C 2 )® 4 ) 
is a subspace of C 2 <g> C 4 (the A subsystem acts on C 2 
and the remaining BCD subsystems act on C 4 ) it holds 
that r(g TA ) < 8. On the other hand, S((C 2 ) m ) is also a 
subspace of C 3 <g> C 3 (the two-party subsystems AB and 
CD act on C 3 ) and hence r(g FAB ) < 9. 

Further constraints on the ranks can be imposed by 
utilizing results of Refs. [l| and [3(J. First of all from 
Ref. [13] we know that any PPT 7V-qubit (N > 3) state 
acting on S((C 2 )® N ) with r(g) < N is fully separable, 
i.e., it can be written as 

e = Ef*e2® ( 41 ) 

i 

This implies immediately that we can restrict our consid- 
erations to r(g) = 5. The same reasoning can be applied 
to g VA and consequently r{g VA ) > 5. Finally, we can 
treat a = g TAB as a three-partite PPT state acting on 
C 2 <g)C 2 ®C 3 (recall that S{<C 2 ®<C 2 ) is isomorphic to C 3 ). 
It follows from Ref. (2j| that if r(g FAB ) < 3 then a is a 
three-partite fully separable (see Eq. (|4"Tjl ) and thus, via 
reasoning from Ref. [30J g is a four-party fully separable 
state. Therefore we can assume that r(g TAB ) > 4. 

Let us now pass to the criterion. The Hcrmitian ma- 
trix we look for has to be also symmetric as we need the 
condition R(H) C R(g) to be satisfied. Consequently 
r(H) < 5 and we have at most 25 parameters to de- 
termine. The condition for r(g TA ) and r(g rAB ) under 
which H exists can be straightforwardly determined as 
in the previous cases and reads [r(g rA )] 2 + [r{g TAB )] 2 > 
121, where one follows from the normalization condition. 
It implies that there are no symmetric PPT extremal 
states with the following triples of ranks (5, 7, 9), (5, 8, 8), 
(5,8,9) (this agrees with the fact that points being in- 
teriors points of two sets and not the extremal point of 
the third one cannot be extremal points of an intersec- 
tion of these three sets). All this implies that we have 
the following ranks for which it is possible that symmet- 
ric four-qubit PPT entangled and extremal states exist: 
r{g) = 5, 5 < r(g TA ) < 7, and 4 < r{g TAB ) < 8. 

Still, however, by virtue of Refs. [l2j and [3(3] further 
simplifications can be inferred. For instance we can treat 
a = g TA as a 2 (g) 6 density matrix with respect to the cut 
B\ACD. Now, if a is supported on C 2 ® C p (3 < p < 6) 
then if follows from Ref. [12| that the condition r(g FA ) > 
p has to be satisfied. Notice that for r{g TA ) = p, a is sep- 
arable across B\ACD and therefore g is fully separable. 
To see this explicitly let us notice that g can be writ- 
ten as g = Z^e® )(ef | \^% b )^% b \ Va . Then, it 
suffices to utilize the fact that tracing out the A subsys- 
tem we get three-qubit symmetric density matrix which 
is separable across the B\CD cut. Finally, utilizing sim- 
ilar approach as the one in Ref. [30] one sees that g has 
to be fully separable. Analogous arguments work also in 
the case of g TAB with respect to e.g. C\ABD cut. 
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In view of what was just said it seems then that 
generically the only possible cases with respect to ranks 
(r{g),r(g rA ),r(g rAB )) are (5,7,7), and (5,7,8). 

VI. CONCLUSION 

Let us discuss shortly the presented results. Once 
more we need to stress that in principle due to results 
of Gurvits @ and Doherty Q the problem of separabil- 
ity in lower dimensional Hilbert spaces can be regarded 
as solved. Still however, it seems that some progress can 
be achieved in systems like 3(g) 3 or 2(g) 4. Here one knows 
that except the states that are detected by the transpo- 
sition map (NPT states) there are also entangled states 
with positive partial transposition. As we do not have 
any unique structural criterion that could detect these 
PPT states, an attempt to learn about the geometry of 
the set of PPT states in the lower dimensional systems 
like 2 ® 4 or 3 <g 3 seems interesting. Some progress in 
this direction in 3 (g 3 systems has recently been obtained 
in Refs. [HI, [H, [I?], H3], where new classes of the edge 
states have been found for all possible configurations of 
ranks of a density matrix and its partial transposition. 
What is even more important examples of these states 
with bi-ranks (4,4), (5,5), and (6,6) were shown to be 
extremal in the set of PPT states in Refs. [S S . 

The main purpose of the paper was to study extremal- 
ity in the convex set of PPT states in qubit-ququart sys- 
tems and the main purpose was to provide a operational 
criterion allowing to judge if a given g is extremal. At 
the proof stage we learned, however, about Ref. [l8| (see 
also Ref. [19( ) in which similar criterion has already been 
given. In this work we have provided our formulation of 
the criterion. It reduces the question of extremality to 
the problem of solving a system of linear equations. Us- 
ing this approach we have reconstructed the known lit- 
erature results, i.e., that there are no extremal states in 
T>2 4 T ranks of which satisfy r(g) + r(g TA ) > 12 (except 
for the case of r(g) = r(g TA ) = 6). We have investi- 
gated the remaining cases of the bi-rank (5,5), (5,6), 



and (6,6). In the case of (6,6) the criterion shows that 
there are no extremal states. In the case of bi-rank (5,5) 
we proved that all the (5, 5) edge states are extremal in 
T>2 4 T implying that the famous Horodecki PPT states 
[5j are extremal. On the other hand, existence of PPT 
entangled extremal states of bi-rank (5, 6) has been con- 
firmed numerically in Ref. [l8l |. however, without giving 
explicit examples. It means that in the 2(g4 systems only 
the PPT entangled states of bi-rank (5, 5) and (5, 6) can 
be extremal. We have also provided a class of states that 
for some parameter region are (5,6) edge states, while 
for some other parameter region they constitute other 
examples of (5,5) extremal states. 

Using the criterion, we have made an attempt to study 
the structure of (5, 6) PPT states. For any (5, 6) g we pro- 
vided such a Hermitian H that R{H) C R(g), while the 
condition R(H Ta ) C R(g TA ) holds if some corresponding 
set of linear equations has a solution. We asked what are 
the conditions imposed on the matrix B (see Eq. (|To| ) 
if the solutions are of some particular form. It is shown 
that whenever B satisfies the corresponding conditions 
the solution to the system of equations exists. 

Finally, we extended our considerations on extremality 
to the simple systems of many-qubits as general three- 
qubit and four-qubit symmetric PPT states. In the first 
case we have shown that the only cases with respect to 
the respective ranks in which one may expect extremal 
PPTES are (4, 4, 4) and (4, 4, 5). 
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